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Abstract. We analyze the Wigner function constructed on the basis of the discrete 
rotation and displacement operators labelled with elements of the underlying finite 
field. We separately discuss the case of odd and even characteristics and analyze the 
algebraic origin of the non uniqueness of the representation of the Wigner function. 
Explicit expressions for the Wigner kernel are given in both cases. 



1. Introduction 

The Wootters' construction of the Wigner function PJ |2j in the discrete phase space 
for quantum systems whose dimension d = p n is a power of a prime number has been 
applied for the analysis of several problems, mainly related to quantum information 
processes IH 13 E] • In this case, the phase space is a d x d grid with a well defined 
geometrical structure (the so-called finite geometry) jjj, which is a direct consequence 
of the underlying finite field (Galois field) structure, so that the phase-space coordinates 
can be chosen as elements of the corresponding finite field PUI2]- The Wigner function 
(a discrete symbol of the density matrix) is constructed using the so-called phase point 
operators and consists in a highly non-univoque procedure of association of states in the 
Hilbert space with some geometrical structures (lines) in the discrete phase space, is a 
substantial disadvantage. 

Slightly different, operational approaches to the Wigner function construction, 
based on the algebraic structure of the Galois fields, have been recently |H1 El EE HI] 
proposed for the quantum systems of prime power dimensions, although, for the prime 
dimensions such constructions were also discussed earlier [Tl| fTTT] . 

In the present article we analyze the Wigner function constructed on the basis of the 
discrete rotation and displacement operators labelled with the elements of the underlying 
finite field. We separately discuss the case of odd and even characteristics and analyze 
the algebraic origin of the non-uniqueness of the representation of the Wigner function 
in both cases. The main difference with the Wootters' method consists in labelling the 
elements of the Hilbert space and the operators acting on this space directly by elements 
of the finite field [12], so that the possible factorizations of the corresponding operators 
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appear as a result of choosing different basis in the field, and are not really necessary 
for the phase space construction. 

The paper is organized as follows. In Sec. 2 we outline the structure of the 
generalized Pauli group and introduce the Wigner mapping as that which satisfies 
the Stratonovich-Weyl postulates and shows the phase problem for the displacement 
operator. In Sec. 3 we analyze the rotation operators for fields of odd and even 
characteristics and discuss the discrete phase space construction. In Sec. 4 we briefly 
outline the reconstruction procedure. In Sec. 5 we discuss different possibilities for 
ordering points in the phase space with elements of the finite field. In Sec. 6 we analyze 
the relation between abstract states and states of physical systems. Several useful 
relations are proved in Appendices. 



2. General definitions 



In the case of prime power dimensions, d = p n , the phase-space representation can be 
constructed in a similar way as in the prime dimensions P^EI]- ^ e P er f° rm the Wigner 
mapping in a slightly different manner than in [2|, by using the generalized position 
and momentum operators introduced in (T2j (see also |H]). In the case of composite 
dimensions instead of natural numbers we have to use elements of the finite field GF(d) 
to label states of the system and operators acting on the corresponding Hilbert space. 
In particular, we will denote as |a) , a G GF (d) an orthonormal basis in the Hilbert 
space of the quantum system, (a |/3) = 5 a ^. Operationally, the elements of the basis 
can be labelled by powers of primitive elements (see Appendix A). These vectors will 
be considered as eigenvectors of the generalized position operator which belong to the 
generalized Pauli group. The generators of this group, usually called generalized position 
and momentum operators, are defined as follows 

Zp\a) =x(a0)\a), X p \a) = \at + (3) , a,f3eGF(d), (1) 
Z p = Z_p, X p = X.p, (2) 

so that 

Z a X p = x {a(3) XpZ a , 

where x{9) is an additive character [Zj 

~2iri 



X{9) = exp 



-tr (9) 
P 



(3) 



and the trace operation, which maps elements of GF (d) into the prime field GF (p) ~ Z p 
is defined as 

tr(0) = 9 + 6 p + e p2 + ... + e pn ~\ 

this operation leaves the elements of the prime field (see Appendix A) invariant. The 
characters (JHJ) satisfy the following important properties: 

X {ap) = d5 ,p, X (« + P) = X («) X W) , (4) 

a^GF(d) 
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similarly to the well known identities in the prime case 
p-i 

u mn = p6 n , , u m+n = u m u n , 

m=0 

where ui = exp[2iii/p\ is a root of unity. 

In particular, for the prime dimensional case (d = p), the position and momentum 
operators act in the standard basis \n),n = 0, ...,p — 1, as follows |TH] . 

Z\n)=u n \n), X\n) = \n + 1), ZX = uXZ, (5) 

where all the algebraic operations are on mod(p), n G Z p , and the rest of the operators, 
analogous to Z a and Xp, can be obtained as powers of Z and X. 

It is worth noting that there is a single element of this basis (sometimes called 
stabilizer state ), labelled with the zero element of the field, which is a common eigenstate 
of all Zp with all the eigenvalues equal to the unity: 

Zp |o> = |o> , 

for any (3 G GF (d). 

The operators (0) are related through the finite Fourier transform operator 



F = \ E x(aP)\a)(P\ t FF* = F*F = I, (6) 

* a,f3£GF(d) 

so that 

FX a F^ = Z a , (7) 

and F A = I for d = p n where p ^ 2, and F 2 = I for d = 2 n . The Fourier transform 
offers us the possibility of introducing the conjugate basis, which is related to the basis 
| a) as follows 

\a) = F | a) , Zp | a) = |a + (3), Xp \a) = x* (a/5) |5) , (8) 

so that the elements of the conjugate basis are eigenvectors of the momentum operators. 

Operators Z a and Xp are particular cases of the so-called displacement operators 
which, in general, have the form 

D(a,(3) = ( j ) (a,(3)Z a Xp, (9) 

and the phase factor <f>(a,/3) is such that the unitary condition 

D{a,P)Di (a,P)=I, 

is satisfied, implying that 

(a, /3)<f («,/?) = 1. (10) 

The operational basis (jOJ becomes orthonormal, 
Tr [D {ai,p 1 )D{a 2 ,fo)] = d6 

— Oil jt*2 



Geometrical approach to the discrete Wigner function 4 

or equivalently D^(a,/3) = D(—a,—(3), if the phase </>(a,/3) satisfies the following 
condition, 

<l>(a l ,P)<P(-a,-P) = x(-al3), (11) 

and for the particular case of the fields of even characteristic, char (GF (d)) = 2, implies 
that 

2 («,/?) = xM), (12) 

which is equivalent to (a, [3) = D (a, (3). 

The displacement operators are non-Hermitian and thus, cannot be used for 
mapping Hermitian operators into real phase-space functions. Nevertheless, a desirable 
Hermitian kernel can be defined as the following transformation [HJ EI] of the 
displacement operator (JHJ), 

A(a,/3) = ~ x(ocX-/3k)D{k,X). (13) 

K,\eGF(d) 

This operator can be used for mapping operators into phase-space functions in a self 

consistent way and satisfies the Stratonovich-Weyl postulates [TB] : 

hemiticity 

A(a : P) = A\a,P), (14) 

if the condition (fTTj) (or (fi"2"j) in the case p = 2) is satisfied; 
normalization 

- d J2 A («^) = / ; 

o,/3£GF(d) 

covariance 

D («, A) A (a, (3) D ] (k, A) = A {a + k, (3 + A) , 
and the orthogonality relation 

Tr (A(a,/3)AV, /?'))= d5 a ,a'W- 

Since the Stratonovich-Weyl postulates are satisfied, the symbol of an operator / 
is defined in the standard way 

W f (a,l3)=Tr[fA(a,f3)}, (15) 

where Tr means the operational trace in the Hilbert space and the inversion relation is 

f = \ E W f (a,P)A(a,(3). 

The overlap relation has the standard form 

Tr(fg) = d W f (a,f3)W g (a,P), 

a,/3£GF(d) 
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and as a particular case, the average value of the operator / is calculated as 

(f) = d W f (a,P)W p (a,P), 

a,/3£GF(d) 

where p is the density matrix. 

In terms of the expansion coefficients of the operator / in the operational basis (J^J) 

/= f«,pD(a,P), (16) 

a,/3eGF(d) 

the symbol of / has the form 

W f (a,P)= Yl /«,aX*M-/^)- ( 17 ) 

K,\eGF(d) 

As some simple examples, we obtain that the symbols of the operators Z K and X\ are 

W Zk {<*,P) = x(P*), W Xx (a,P) = X (-<xX), (18) 
and in the particular case of prime dimensions we get for the symbols of (jSJ) 

W z (a, b) = u\ W x (a, b) = u~ a , 
where a, b £ Z p . In the same way, the symbols of the basis states \k) and \k) are 

W\ K )( K \ (a, P) = 5p tK , W\~ k){ ~ k \ (a, (3) = 5 a>K . 

3. Discrete phase space construction 

3.1. Lines and rays 

In the discrete space GF (d) x GF (d) the concept of line can be introduced in a similar 
way as in the continuous plane case, so all the points (a, 0) £ GF (d) x GF (d) which 
satisfy the relation 

(a + V (3 = ti, 

where (, rj, 9 are some fixed elements of GF (d), form a line. Moreover, two lines 

£a + 7ip = ti, ('a + ri'P = (19) 

are called parallel if they have no common points which implies that 77C' = C 7 !' '■ If th e 
lines (|T9*|) are not parallel they cross each other at a single point with coordinates 

a = w - r/tf) (Cv - CvT 1 , P = (C^ - C'tf) iCv' - Cv)' 1 ■ 

A line which passes through the origin is called ray and its equation has the form 

a = 0, or p = fiat (20) 

so that a = and P — are the vertical and horizontal axes correspondingly. 

Each ray is characterized by the value of the "slope" \i and we denote by A M the 
ray which is a collection of points satisfying P = (ia and by Aoo the ray corresponding 
to the vertical axis. 

There are d—1 parallel lines to each of d + 1 rays, so that the total number of lines 
is d(d + 1). The collection of d parallel lines is called striation j2]. 
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3.2. Displacement operators 

The displacement operators labelled with points of the phase space belonging to the 
same ray commute (here we omit the phase factor): 



X ( (1>CXiCX2) Z ai + a2 X ^( ai +a 2 ) Z a2 X f) 2= ^ a2 Z ai X p 1= ^ ai , 

and thus, have a common system of eigenvectors {|V#)> I 1 , v e (^)} : 

Z a X^|<) =exp(^)|<), (21) 

where /x is fixed and exp(i^ /ljl/ ) is the corresponding eigenvalue, so that j^) = are 
eigenstates of Z a operators (displacement operators labelled with the points of the ray 
(3 — - horizontal axis) and = F = |z>) are the eigenstates of Xp operators 
(displacement operators labelled with the points of the ray a — - vertical axis). 

In the simplest cases, d = 3 and d — 4, the displacement operators have the form 
(up to a phase) 

ray equation displacement operators ray equation displacement operators 
P = Z,Z 2 (3 = Z e ,Z e 2,Z d 3 

(3 = a ZX, Z 2 X 2 (3 = a ZqXq, Z^Xqz, ZqzXqz 

13 = 2a ZX 2 ,Z 2 X 4 = Z 2 X (3 = 6a Z e X 6 2, Z 2X e z, Z e zX e 
a = X,X 2 (3 = 6 2 a Z e X ffA , Z 9 2X , Z ffA X ffA 

a = Xg, Xg2,X e :i 

where is the primitive element, a root of the polynomial 6 2 + 9 + 1 = 0. 

It is easy to observe that an arbitrary displacement operator Z T X V acting on 
an eigenstate of the set {Z a X^ a , jj, is fixed, a G GF (d)} transforms it into another 
eigenstate of the same set: 

Z a X m [Z T X V \^)\ = exp(i^ u )x(av - fiar)Z T X v ^). 

It is clear that if v = fir we do not generate another state (in other words we do 
not change the index is). Because for an arbitrary v one can always find such x that 
v = fir + x, we can generate all the states from the set {\ip£),(J> is fixed} by applying only 
the momentum operators X v (where v runs through the whole field) to any particular 
state belonging to this set. It is clear that all the eigenstates of Xp operator can be 
obtained by applying Z a operator to any particular state from the set { 



€ >}. 



3.3. Rotation operators 



The "rotation" operators which transform eigenstates of the operators associated to 
the ray (3 = /ia: 

{/, Z ai X^ ai , Z a2 X^ a2 , ...} (22) 
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into eigenstates of the operators labelled with points of the ray /? = (// + ji') a: 

{I, Z ai X^ + ^) ai , Z a , 2 X^ + ^) a2 , ...} , (23) 
are defined through the relations 

VpZ a Vt = exp(vp (a, ii))Z a X m , [V^ X u ] = 0, V = I, (24) 

for all fi, v e GF (d). 

Really, being (|2*T|) a state assigned to the ray f3 = fia, we obtain after simple 
algebra 

v^z a x^} = v^z a x m (y}fVj) K) = 

exp(^)V^K) = exp(iip(a, /i'))Z a X (f , +lxl)a V^\^) } 

that is 

Z a X^ )a \VM)] = exp(i(^,, - <P («,//))) t^'IO] > (25) 

i.e. the state V^ipy) is an eigenstate of the set (J2HJ) - This means that we can interpret 
the action of operator as a "rotation" in the discrete phase space: 

-V ~~* V+m'> (26) 

although it should be careful in the case of fields GF(2 n ) as we will discuss below. 
Observe, that one cannot reach the vertical axis by applying to any other ray. 

The explicit form of can be found taking into account that they are diagonal in 
the conjugate basis (|2*4*|): 

V^= 22 c «,mI^)(^I' c o,m = !• (27) 
Transforming the position operator Z a with V^, 

VpZ a V*= c ^+^ c lJ K + a )(^\ 

K&GF{d) 

and taking into account that 

Z a X fia = X (-fictK) \k + ot)(«|, 

KdGF(d) 

we obtain that the coefficients c K satisfy the following condition 

c K+a^c* K ^ = exp(iip (a, fj))x (-/ia«) . 
In particular, for k = we obtain 

exp(iy2 (a, //)) = c ajM c^ = c Q)M , (28) 

that is 

c K+ a ^c* = c a ^x , (29) 
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and substituting a = we get |c Kj/ J 2 = 1, which also follows from the unitary condition 

vrf = v}v, = I. 

It is easy to note that Eq. (j2HJ) is automatically satisfied after the substitution 

c a,fi * C a,fi c a,^X ( Ctv) , 

which means that all different sets of operators have the form 

= V^X U , (30) 
where is constructed using an arbitrary solution of Eq. (|2T?j) . 

3.3.1. Fields of odd characteristic In the case of fields of odd characteristic, we impose 
an additional restriction on the rotation operators: we demand that V M form an Abelian 
group: 

= w> ( 31 ) 

which implies that c K ^c K ^> = c KjAt+M / and in particular c* jAt = c Kj _ M , leading to the 
relation V£ = V^^. In this case the relation (}2l)j) is well defined, i.e. the operator 
transforms a state associated with the ray A^ into a state associated with the ray A M+At '. 
It can be shown (see below) that the condition (|31j) cannot be satisfied for the fields of 
even characteristic, so that this case should be considered separately. 
Then, the solution of Eq. (|29|) can be easily found 

c K>ft = X (-2~ V/i) , (32) 

so that 

V » = £ x(-2 _1 «V)l«>(«l- (33) 

K£GF(d) 

In the prime field case, GF(p), p ^ 2, the whole set of rotation operators is produced 
by taking powers of a single operator [THj 

p-i 

V = J2^{-2~ 1 k 2 ) \k)(k\. 

k=0 

In particular, the state V™']^™) is associated with the ray b = (m + m') a, where the 
algebraic operations are modp, so that 

Ao — >• Ai, Aq — > Ai — > X2 , etc. 

v 2 

3.3.2. Fields of even characteristic The situation is more complicated for fields of 
char(GF (d)) = 2. Really, it follows from (|27IJ) that (substituting k = a) 

= X • (34) 

The solution of the above equation is not unique and thus, there is an ambiguity in 
solving Eq. (f2~H|) (see Appendix B). 
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One of the consequences of this ambiguity is that the operators of the form (|27|) . 
where c Kj/J is a particular (for a fixed value of fi and k G GF(2 n )) solution of (|2"§j) do 
not form a group. In particular, the operator V 2 is not the identity operator. Really, 
using we have 

K€GF(2 n ) K£GF(d) 

which due to the property tra = tra 2 , a G GF (2 n ), can be transformed into 

V l= E X^/i 2 "" 1 ) |J5)^l=^i, (35) 

kGGF(2") 

where the relation n 2 fi = ^k/z 2 ™ j has been used. 

It also follows from (jHEJ) that inside the set {V^,/x G GF (2 n )} an inverse operator 
to a given from this set does not exist. To find the inverse operator to some we 
have to extend the set {V^, \i G GF (2 n )} to the whole collection of all the possible 
rotation operators defined in (jHHJ), i.e. to the set {V^ u , /J, v G GF(2 n )}. Then, it is 
easy to conclude from (|33j) that 

which implies the following relation between c a ^\ 

2 n-l 



c*a,n = X[a^ J Co, M . (36) 

We will fix the operator V^ ;I , = o is such a way (see Appendix B) that the coefficients 
c K)At , corresponding to the basis elements of the field, k = a\, a n , are chosen positive, 
so that we have 

(n— 1 n \ 

nJ2 ki(7i E k 3 a 3 U ?=i\/x(kfaffi), (37) 
i=l j=i+l / 

n 

k = y^^kjCTi, h G Z 2 , (38) 
i=i 

where {<Tj, i = 1, ..n} are the elements of the basis and the principal branch of the 
square root in (j3*7|) is chosen. All the other possible rotating operators can be obtained 
according to (QUI) . 

The whole set of operators {V^ u , //, v G GF (2 n )} form a group. Really, let us pick 
two operators and constructed according to (pT7j) . Using the properties and 
(|2"%j) we obtain 

On the other hand we have 

which suggests that (recall that |c QiM | = 1) 

c a ^c a ^ = exp (if (a, /i, //)) c a ^ +fl >, (40) 
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where f(a,fi,fi') = f(a, //, /i) is a real function of and //. Let us note that due 
to the property (j2HJ), the function /(a,/i, //) linearly depends on the parameter a, in 
the sense that f(a + f3,fi,fi') = f(a,fi, //) + /(/?,//,//). The complex conjugate of 
(HOJ) together with (pjj) leads to the condition exp (if (a, /i, //)) = exp (— i/(a, /i, /i')) = 
±1 G ^2- Because every linear map from a G GF(p n ) to GF(p) ~ 2p is of the form 
of a character 7j: a — > x( a P)^P G GF(p n ), the function exp (if (a, /i, //)) can be 
represented as 

exp (i/(a, fi'))=x (af(ji, //)) • 

It follows from (|40|) and from the above equation that there exists z/ = /(/i, //) (see 
Appendix C) so that 

V,y„Z a Vffl = X v V^Z a Vl^Xl = c a ^x K) Z a X {fl+IM , )cn (41) 

which immediately leads to the relation 

Finally, in general we have 
3.4- Phase space construction 

We can associate the lines in the discrete phase space with states in the Hilbert space 
according to the following construction: 

1. The eigenstate of {Z a } operators with all eigenvalues equal to 1, |0) = 1^°), 
(note that such state is unique) is associated with the horizontal axis, (3 = 0. It is worth 
noting that such an association is arbitrary and in some sense fixes a definite quantum 
net 0. 

2. All the other states of the "first" striation are obtained by applying the 
displacement operator X v to |0), so that the state \ip®) = X u \0) is associated with 
the horizontal line which crosses the vertical axis at the point (0, z/), i.e. with the line 
[3 = i'. The states are eigenstates of the set {Z a }: 

Z<M) = x(<*v)X v Z a \0) = x(<xuM), (42) 

and form an orthonormal basis (V^IV^') = <W. 

3. All the other striations are constructed as follows: First we apply the rotation 
operator to the state |0) and the obtained state j^o) = V^|0) is associated to the 
ray (3 = \ia. The state [iptf) is an eigenstate of the set {Z a X^ a } according to (J25J) . It 
is worth noting that different sets of rotation operators can be chosen, which leads to 
different associations between states and lines (see discussion below). 

4. All the other states of the /z-th striation are obtained by applying the operator 
X u to the state \4>q): 

\W) = X V \1$). (43) 
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The states \ip„) are eigenstates of the set {Z a X^ a }: 

ZaX^ltf) = expHV (a, fx))V^Z Q V^X u V^\0) = 

exp(-iip (a, fi))V^Z a X u \0) = x(az/) exp(-i(p (a, ^))X V V^}) 

= X(ajy) exp(-i(p 

and are associated to the lines (3 = \ia + v. The phase exp(— i<p (a, f/,)) in the above 
equation is defined in pij). 

Note that \ip^) = V^ u \0), where V^ u is defined in flSUj) . As it was shown in [T2j the 
states (jHH) form mutually unbiased bases (MUBs) jT2j: 

l«l€')l 2 = ^ m^m', 

which have been extensively discussed from different points of view in recent papers 
[TH1 Uni H2 UJ] • In particular, it is known [3 E| that different MUBs can be associated 
with different striations. 

In our approach we note that the inner product of two states associated to lines 
can be rewritten as follows 

(WW') = {OtyXp^VM = (0\X^V^\0) = J2 c k ,^{v-i/\k){k\0), (44) 

K&GF{d) 

so that, taking into account = x( Kh> ) I 'd 1 ' 2 '', we obtain 

KdGF(d) 

Then, we get 

i*s£r= E ^'-m^-mx((«-«o(^-o), 

K,n'£GF{d) 

so that after changing the index k — k! = A, and making use the relation (J29j) we 
immediately obtain 



K',\eGF(d) 




The uppermost relation corresponds to the inner product of states belonging to the 
same striation, i.e. to the same basis, and from the other one can observe that different 
striations correspond to different MUBs. 

5. The state associated with the vertical axis (3 is obtained from |0) by applying 
the Fourier transform operator: 

|0)=F|0). (45) 

The state (|4*5)l is an eigenstate of the set {Xp, [3 G GF(d)} with all eigenvalues equal 
to unity: X^O) = |0). The vertical lines, parallel to the axis (3, crossing the axis a at 
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the points a = u are associated with the states \ip®) = Z u \0), which are eigenstates of 
the set {Xp}: 

XpHl) = X (-(3v)ZvXp\0) = X {-^M)- 
It is clear that 

WM)\ 2 = \- 

For the fields of odd characteristics we can always choose the set of rotation 
operators in such a way that they form a group, which in particular means that a 
consecutive application of any from this set returns us to the initial state: V£ = I. 
This does not hold for the fields of even characteristic due to Eq.(|35|). which can be 
interpreted as follows: the first application of the rotation operator transforms a 
state associated with the ray of some striation to a state corresponding to a ray from 
another striation. The second application of does not return to the initial state (as it 
could be expected from (|26)l). but to a state which belongs to the original striation and 
displaced by X 2 n-i. So that, from the geometrical point of view the operator for 
GF(2 n ) transforms a ray from some striation to a parallel line from the same striation. 

3.5. The phase of the displacement operator 

The phase <p (a,p) which appears in (J24*j) is intimately related to the phase <fi (a>,f3) of 
the displacement operator (jSJ). Let us impose the natural condition that the sum of the 
Wigner function along the line (3 = pa + v is equal to the average value of the density 
matrix over the state associated with that line £Q: 

- d Y, W(a,P)6p, m+u = {W\p\W), (46) 

a,/3eGF(d) 

where 

W(a,P) = Tr\pA(a,P)]. (47) 

According to the general construction, the state is obtained by application of the 
operator V^ u to the state |0): = V^ u |0). After simple algebra we transform the 

left-hand side of (|4*d^) as follows 

- W (a, P) Sfrpa+v = 2 Yl W(a,pa + v) (48) 

a,(3eGF(d) a£GF(d) 

= 2 2 P«,A0(^ 1 (-^ + A),-/t + A)x(/i~ 1 (^ + A)(A + z/)), 

K,\aGF(d) 

where p = ^ K a Pk,x\^) (A|, meanwhile the right-hand part is converted into 

mpm = (o\v}xipx v v,\o) 

= J? Pr,vC* K ^C Kf ^x(-^(r-u)+K'(v-u)). 
k,k',t,v 
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Taking into account the relation (|29|). the above equation can be rewritten as follows 
(€1 P K> = \Y1 ^VH-^ 0" _1 (« + A) (A + !/)) . (49) 

re, A 

Comparing PH|) and (J4T?|) we observe that (/i -1 (— k + A) , — k + A) = c M -i(_ re+A w, or 
in a compact form 

(r,i>) = c T>T -i v . (50) 

Also, we impose the conditions (r, 0) = 0(0, i>) = 1, which means that the 
displacements along the axes a and 0, which are performed by applying Z K and Xa 
operators correspondingly, do not generate any phase. 

This immediately leads that for fields of odd characteristic 

<P(t,v)= X {-'2- 1 tv), (51) 

and the condition (fTTj) is automatically satisfied. Thus, the displacement operator in 
this case has the form 

D{a,0) = X (-2- 1 aP)Z a X f3 (52) 

so that the unitary condition (a, 0) = D (—a, —0) is satisfied and the kernel (fH?j) 
can be represented in the familiar form 

A(a,0) = D(a,l3)PD ] (a, (3), 

where 

P = - d D (®>P), P\a) = \-a), 

a,f3£GF(d) 

is the parity operator. 

In general, for a density matrix defined in the standard basis \ct) as 

P = \P) (H , (53) 

/j,,uGGF(d) 

the Wigner function takes the form 

Wp ^' ® = d ^ X ^ ^ ~ ® + a ^ ~ ^ ^ ~ p w ^ 

which still can be simplified (summed over 7) for the fields of odd characteristics taking 
into account the explicit expression ()51|) for the phase 0. Nevertheless, for the fields 
of even characteristic, d = 2 n , there is a freedom in the election of the phase 0, which 
takes values ±l,±z, related to different possibilities of choosing the rotation operators 
Vfi,u = VfiX u for the phase space construction. Once fixed the set {V^} of rotation 
operators we can find the phase from (|5()jl and thus, construct the Wigner function 
(El. 
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3.6. Non-uniqueness of the Wigner functions 

It is worth noting that the Wigner function constructed using different distributions 
of signs in V M operators are not trivially related to each other. This difference is of 
fundamental significance for fields of even characteristic (because there is no natural 
choice of the set of rotation operators), although similar considerations can also be 
taken into account for fields of odd characteristic (if different to (p£Ij) set of rotation 
operators is used for the phase-space construction). 

In the rest of this section we will focus on fields of even characteristic. Let us fix 
operators according to (|37|) and choose some other set, so that V^hfjt) — V(iXh(u)> 
where h(n) is an arbitrary function satisfying the conditions: a) h(0) = 0, which basically 
means that even after changing the rotation operators, the displacement operators along 
the axes a and (3 do not have phases, D (a, 0) = Z a , D(0,(3) = Xp; b) the non- 
singularity: ah (a -1 /?) = 0, if a = for any (3 G GF(2 n ), which implies that the 
coefficient co )At in (}27j) is fixed, co >At = 1, for all V^. The simplest example of such 
function can be given in the case when rotation operators are labelled with powers of a 
primitive element a: {V^fJ, E GF (2™)} = {V a k,k = l,..,2 n - 1}, then h(a k ) = a n < k \ 
where m is a natural number which depends on the value of k. 

The Wigner function constructed using the new rotation operators V^ v has the 

form 

W' p (a, 0) = Tr [pA' (a, (3)] , A'(a,/3) = 1 £ x (aA - (3k) D' («, A) , 

where 

D' («, A) = c' K K - lx Z K X x , 
and d K £ are the matrix elements (|2*7jl of V p ^ p ) in the conjugate basis \a), so that 

C U-!A = X (Kh (k _1 A)) C k>k -i\. 

Then, we have 

D' (k, A) = x { K h («T A)) c KyK -i\Z K X x = c K)K -\\Xh( K -i\)Z K Xh( K -i\)X\ 
= X h ( K -ix)D («, A) X h ( K -ix). 

After some simple algebra we obtain the "new" Wigner function W' p (a, (3) in terms of 
coefficients of the expansion 

W ' P ( a , P) = E p ^ uX ( a ( u ~i I )-fl7 + l h {t 1 iy ~ /")) + iv) <f>[rf,v- 

which is related to the "old" Wigner function (constructed with V p operators) as follows 
W' p (a, (3)=^ J2 W p (a + 1 h (7- 1 ) + 77, P + i) (56) 

7' 6GF(2") 

+ ^ E W p( a >P) + ^ E (57) 

/3GGF(2«) <*eGF(2") 
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Note, that 

/3eGF(2") aeGF(2 n ) 

are the marginal probabilities to detect the system at the states |a) and \$) 
correspondingly. 

It is clear that the sum of the "new" Wigner function over a line (3 = pa + v gives 
the same result as the sum of the "old" Wigner function over the points of the shifted 
line P = pa + v + h (p): 

^2 W(a, p) 5p =im+u = ^ W(a,(3) 5p =m+v+h{ ^ . 

a,(3eGF(2 n ) a,/3eGF(2 n ) 

For a particular choice h (p) = k - const, h(0) = 0, we find that the new Wigner 
function has the form 

W' p (a, 0) = W p (a, P + k)- Pp+ K ,0+ K + p p ,p. (58) 

For a more complicated case h (£) = k£ the Wigner function acquires the following form 

W' p (a, P) = W p {a + k, p) + ^p(a) - ^p(a + k). 

Using the relation (J47|) we can find the Wigner function for the state (|43|) 
corresponding to the line P = pa + v. First of all, using (|27|) and the property 
(ac|A) = x(— K\)/d 1 / 2 we rewrite the state as follows 

^ = ~n Yl c wx(- KV )\*) = 2 C ^M K ( X ~ z/ ))l A )- 

K£GF(d) n,\£GF(d) 

Then, from (|T5jl and (JTTjl we obtain after long but straightforward algebra 
W m (a,P) = «|A(a,/5)|C) = <W„, 

for fields of both odd and even characteristic independently on the choice of the rotation 
operator in the last case. 

In general, the Wigner function of any state of the form (}53|) with p^ u = q^5^ u 
does not depend on the sign choice (and thus, is uniquely defined) as it can be observed 
directly from (jB3j) . 

As a nontrivial example of essentially different Wigner functions which can be 
associated to the same state let us consider a particular case of pure states with real 
coefficients, i.e. p^ u = p v ^. For this class of states the Wigner function (jB^|) can be 
rewritten in the following explicitly symmetric form 

K = ^PT PwX(a(v + p)+Pl + lh(j- 1 (is + p))) 

77^0 

x0(7,^ + /^) IxM +X(1P)} + 7^ Yl PwX(a(v- p)) +P/3./3- 

^udGF(2 n ) 
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Now, we observe from the above equation that if the factor [x (jf) + X (if 1 )} is zero, 
the Wigner function obviously does not depend on the choice of the function h{y). 
Thus, two Wigner functions corresponding to functions h\{v) and hi iy) are the same 
if simultaneously x (lh (7- 1 (u + //))) = x (7^2 (l~ l (f + fi))) and x (7^) + X (in) ^ 
or, in other words, 

tr (7 [hi (7- 1 (u + 11)) + h 2 (7- 1 {u + fi))]) = 0, (59) 

tr (7 [// + *,]) =0, (60) 

where 7 / and v. 

For instance, consider the state \ip) = (|0) + |cr 3 )) /y2 in the case GF (2 2 ) and fix 
the irreducible polynomial as cr 2 + cr + 1 = 0. The indices \i and 1/ in (|59* |) -(|5l1 J) take 
on values and a 3 , so that 11 + v = a 3 , and thus, it follows from (J60|) that the only 
admissible value of 7 is 7 = cr 3 . Then, the condition (JHHJ) leads to the following equation 



for the functions h\{y) and h 2 (u): 

hi {a 3 ) + h 2 (a 3 ) = a 3 . (61) 

The first set of solutions is h\ (a 3 ) = a, hi (a 3 ) = cr 2 and correspondingly hi (a 3 ) = cr 2 , 
hi (cr 3 ) = cr. This means that two sets of rotation operators 

1. X K V a ,X x V a 2,X a V a3 , (62) 

2. X^XyV^X^, (63) 



where k, k', A, A' G GF (2 2 ) lead to the same Wigner function. Note, that the constants 
K, k', A, A' are arbitrary elements of GF (2 2 ) because there are no restrictions imposed 
neither on h (cr) nor h (cr 2 ). 

The second set of solutions of (jBTj) is hi (a 3 ) = a 3 , h 2 (cr 3 ) = and, correspondingly, 
hi (cr 3 ) = 0, hi (cr 3 ) = cr 3 , so that the rotation operators 

3. X^XxnV^X^V^ (64) 

4. X^V^XxnV^Vo*, (65) 

with k", k'", A", A'" G GF (2 2 ) produce the same Wigner functions. 

Finally, there are only two different Wigner functions to represent the state 
1^) = (|0) + I " 3 )) /\/2 (compare with 0). It is worth noting that this state is labeled 
by the elements of the field GF (2 2 ) and thus has no direct relation to the physical state, 
until the basis for the field representation is fixed (see Sec. 6). 

Obviously, for larger fields the variety of Wigner functions for a given state rapidly 
grows with the dimension of the field even for highly symmetrical states. For instance, 
the state = (|0) + |cr 7 )) /y/2 labelled with elements of GF (2 3 ) can be represented in 
8 different ways. 

For fields of odd characteristic, in the whole group {V^X U , //, v G GF(p n )} a 
subgroup containing only rotation operators {V^, 11 G GF(p n )} can be separated, which 
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allows us to construct the phase space as it is outlined in previous sections, so that the 
Wigner function is uniquely defined for a given state. Nevertheless, the whole group can 
be used for phase-space construction as well, which would lead to the non-uniqueness 
in the definition of the Wigner function, very similar to ([56)1 . 

For an arbitrary state, we can easily calculate a total number of possible Wigner 
functions which represent this state in the discrete phase space. According to the present 
construction we fix the phase of the state corresponding to the horizontal line ([42)1 . 
Besides, we fix the property (jJJ) of the Fourier transform operator (JHJ), i.e. the Fourier 
transformation of Z a operators generate X a operators without any phase factor (which 
in principle is not necessary if the property F A = I for d = p n , where p ^ 2 and 
F 2 = I for d = 2 n is not required). Now, we can generate all the possible Wigner 
functions choosing different (d—1) rotation operators V^X V (both for fields of odd and 
even characteristics), which gives d d ~ x different structures (which is directly related to 
different quantum nets introduced in 0). Nevertheless, the symmetry of the state can 
essentially reduce the number of different Wigner functions. 

4. Reconstruction procedure 

It is well known that the Wigner function can be reconstructed using projective 
measurements, associated with a summation over the lines PQI2DI- in this section we 
explicitly relate the elements of the density matrix with the corresponding tomogram, 
the averages of the form (if)$\p\if)$) = u;(/x, v). 

Let us start with the relation of the tomogram u>(p, v) with the Wigner function: 

a,PeGF(d) a£GF{d) 

and consequently with the components of the density matrix in the basis of the 
displacement operators (fTHjl. which is obtained by taking into account the relation (fTTjl : 

K£GF(d) 

The above equation can be immediately inverted using (jlj): 

ueGF(d) 

or, changing the indices 

v£GF(d) 

where 

To reconstruct the matrix element po,A we have to use the results of measurements 
in the conjugate basis \k): 

u (k) = (k\p\k) = 2 Y w ( a ,P)^,K= Y Po,ax(-kA), 

a,/3£GF(d) XeGF(d) 
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which leads to 

On A = 

d 



5. Ordering the points on the axes 

In order to plot the Wigner function we have to choose an arrangement of the field 
elements and take it into account to fix the order of the elements on the axes in the 
finite plane GF (d) x GF(d). When the dimension of the system is a prime number 
there is a natural ordering of the elements, i.e. for any a, b G Z p , a ^ b there is a relation 
either a < b or a > b, so that the points on the axes can be enumerated. For instance, 
in the case Z 5 the elements are arranged as {0, 1, 2, 3, 4}. 

For field extensions there is no natural ordering of elements, however, there are 
several possibilities to arrange elements of the field. The most common choice is an 
ordering according to powers (which are natural numbers and thus can be naturally 
ordered) of some primitive element. Nevertheless, the primitive element, in general, 
is not unique. Actually, for the d- dimensional system there are <fi(d— 1) primitive 
elements, where <fi (r) is the Euler's function, which indicates the number of integers 
with 1 < n < r which are relatively prime to r [7j. 

As an another possibility we can enumerate the points on the axes according to 
the following procedure: first, we fix that the origin corresponds to the point (0,0); 
second, we choose some basis {o~i, o~ 2 , ■ ■ ■ , o~n} in the field GF(p n ) and expand elements 
in this basis: a = YTj=i a j a ji a j e Now, we arrange the expansion coefficients in a 
sequence to form a number on the base p (binary, ternary system, etc.) starting with 
the coefficient a n and taking the leftmost place, it will be followed by a n _i and so on; 
obviously, ai takes the rightmost place. The full number is (a„a n _i . . . a\) p and can 
be transformed into some integer in the decimal basis using the standard procedure: 
(a\ x p°) + (a 2 x p 1 ) + . . . + (a n x p n ~ l ). This procedure is obviously not unique due 
to the existence of different bases in the field. As an example, let us take GF (2 3 ), the 
choice of the irreducible polynomial as x 3 + x + 1 = 0, and use the normal self-dual basis 
{a 3 , a 6 , cr 5 }, in this basis the field elements are 

e = e 6 + e 5 , e 2 = e 3 + e 5 , 

e 3 = e 3 , e 4 = e 3 + e\ 
e 5 = e 5 , e 6 = e 6 , 
e 7 = e 3 + e 6 + e 5 , o. 

Now, we can associate each element with a number in the binary system 
- (110) 2 - (6) M , £ 2 ^(101) 2 ^(5) 10 , 
03 - (001) 2 - (l) u , 4 ^(O11) 2 ^(3) 1O , 

05 _> (100)a _> (4)iQ j £6_ (010)2 ^ (2)io) 

^(111) 2 ^(7) 10 , ^(000) 2 ^(0) 10 . 
So, we arrange the points on the axes using the above ordering according to 

{o,e 3 ,e G ,e A ,e 5 ,e 2 ,6,6 7 }. 
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Another possibility to arrange the elements of the field is in accordance with the 
value of the trace of each element (which is a natural number) ^U] and inside the set 
of the elements with the same trace we can use any of above mentioned ordering, say 
powers of a primitive element. 



6. From abstract states to physical states 

In applications we have to establish a relation between abstract states labelled with 
elements of the field and states of a given physical system. Such interrelation strongly 
depends on the character of a system, for instance, if the system is actually a single 
"particle" with p n energy levels or it consists in n "particles" (degrees of freedom) with 
p energy levels. In the last case the mapping Hd TC P (g> H p .. <E> 7~C P from the abstract 
Hilbert space to n-particle vector space can be achieved by expanding an element of the 
field in a convenient basis {a±, . . . , a n }\ a = a\o~i + . . . + a n a n , a, G Z p , so that 

|a> -> \a 1 } 1 ® . . . ® \a n ) n = \a u . . . , a n ) , 

and the coefficients aj play the role of quantum numbers of each particle. For 
instance, in the case GF{2 2 ) the state (|0) + |cr 3 )) /y/2 corresponds to the physical state 
( 1 00) + 1 10} ) /y/2 in the polynomial basis (l,er), whereas in the self-dual basis (o~,o~ 2 ) it 
is associated with ( 1 00) + |11)) /y/2. Observe, that while one state is factorizable, the 
other one is entangled. 

This implies that all the operators Zp are factorized into a product of single particle 
Z operators, Zp = Z bl ® . . . eg) Z bn , where (3 = b' x o' x + . . . + b' n a' n and {a[, . . . , a' n } is 
the basis which is dual to {o~\, . . . , a n }, see Appendix A, and b\ e Z p . To have a better 
understanding of this aspect let us recall the definition of the operator Z a 

f3eGF(d) 

and choose a basis . . . , a n } to expand (3. Then, taking a = a[ as an element of the 
dual basis, we obtain 

p- 1 / 27r ; \ 

z < = E x W) (A = n u E ex p ( — h ^ (^i) ) • 

Now, if a) i ^ j the duality means tr (c^CTj) = 0, and thus 

E exp ( 6,-tr (ofo) J |6j-) = E \ h i) = J i' 

where the index j means the j-th particle; 
b) i — j we have tr (cr^Cj) = 1 and then, 

E ex P ( hi! (cr-O-i) J |6i) (6i| = E eX P ( ^ ) ( b i\ = Z i- 

bi=0 ^ P ' bi=0 ^ P ' 
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Finally, we obtain the factorization 

Z a <. = I x <g> . . . <g> Ii_ x ® Zi® . . . I n , 

i.e. in the i-th place is located Z operator and all the others are unities. 

In the case of a single "particle" the states of a physical systems can be labelled by 
the elements of the field arranged in some order (0, ai, a p n-i) (see previous Section). 
Then, the free Hamiltonian of the systems takes on the form 

H = Eq |0) (0| + E\ \a.\) (a.\\ + . . . + E p n\a pn -i) (a p n-i|, 

where the energies are arranged in the non- decreasing order: Eq < E\ < . . . < E pn . 

7. Conclusions 

In this paper we studied an explicit form of the kernel operator (the phase point operator 
PU which maps states of a quantum system of dimension d = p n into a Wigner 
function in a discrete phase space. The crucial point in the phase space construction 
play the rotation and displacement operators labelled with elements of GF(d). These 
operators are resulted explicitly related after imposing the condition (pjj) and allow us 
to establish a clear correspondence between states in the Hilbert space of the system and 
lines in the discrete phase space. The structure of the rotation operators is quite different 
for fields of odd and even characteristic. While for the fields of odd characteristic 
the rotation operators form a p n dimensional Abelian group, the corresponding group 
in the case of GF(2 n ) is of order 2 2n and includes both "rotation" and "vertical" 
displacement operators X^fi G GF{2 n ). So that although, for a particular phase- 
space construction the group property is not really unnecessary, different choices of sets 
of rotation operators, lead to different Wigner functions, which is directly connected to 
the freedom in the election of quantum nets in the Wootters' construction [2]. Such 
freedom obviously exists also in the case of odd characteristics, which nevertheless can 
be avoided fixing the rotation group in a natural way ((221) • 
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8. Appendix A: Finite fields 

A set £ is a commutative ring if two binary operations: addition and multiplication 
(both commutative and associative) are defined. 

A field F is a commutative ring with division, i.e. for any a G F there exists 
a -1 G F so that a~ 1 a = aa^ 1 = I (excluding the zero element). Elements of a field form 
groups with respect to addition F and multiplication F* = F — {0}. 

The characteristic of a finite field is the smallest integer p, so that p ■ 1 = 
1 + 1 + .. + 1 = and is always a prime number. Any finite field contains a prime 

p times 
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subfield Z p and has p n elements, where n is a natural number. Moreover, the finite field 
containing p n elements is unique and is usually called Galois field, GF(p n ). GF(p n ) is 
an extension of degree n of Z p , i.e. elements of GF(p n ) can be obtained with Z p and 
all the roots of an irreducible polynomial (That is, one which cannot be factored in Z p ) 
with coefficients inside Z p . 

The multiplicative group of GF(p n ) : GF(p n )* = GF(p n ) - {0} is cyclic 9 P " = 9, 
9 e GF(p n ). The generators of this group are called primitive elements of the field. 

A primitive element of GF{p n ) is a root of an irreducible polynomial of degree n 
over Z p . This polynomial is called a minimal polynomial. 

The trace operation 

tr(a) = a + a 2 + ... + a p " 1 , 
maps any field element into an element of the prime field, tr : GF(p n ) — > Z p , and 

a tr(a) 

satisfies the property 

tr(«i + a 2 ) = tr(ai) + tr(a 2 ). (67) 
The additive characters are defined as 

X(a) = exp 
and possess two important properties: 

x(ai + «2) = xMxM 

and 

E x(«) = o. 

aeGF(p") 

Any finite field GF(p n ) can be considered as an n-dimensional linear vector space 
and there is a basis {o~j,j = l,..,n} in this vector space, so that any a G GF(p n ), 
a = YTj=i a j a j an d a j ^ ^p- There are several bases, one of them is the polynomial basis 
{1, 9, 9 2 , 9 n ~ 1 }, where 9 is a primitive element of GF(p n ), another one is the normal 
basts {6,0P,...,ffl ' }, so one can choose whatever according to the specific problem. 

Two bases {a±, . . . , «„} and . . . , f3 n } in the same field are dual if tr («j/3j) = 
Sij.A basis which is dual to itself is called self-dual basis, tr (ctjCtj) = 

Example: GF(2 2 ), the primitive polynomial is x 2 + x + 1 = 0, it has the roots 
9 2 }. The polynomial basis is {1,0}, whose dual basis is {9 2 , 1}: 

tr(10 2 ) = l, tr(ll)=0, 
tr(99 2 ) = 0, tr(01) = l. 

The normal basis {0, 2 } is self-dual: 

tr(00) = l, tr(00 2 ) = O, 
tr(0 2 0) = O, tr(0 2 2 ) = l. 



2ni 

tr (a) 

V 
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9. Appendix B: Solution of the equation c re+Q!)/Lt c* = c a ^x (^ock) in the d = 2 n 
case 

To solve Eq. (|2U|) we first fix a basis in the field GF{2 n ): {o~j, j = 1, .., n}, so that any 
element of the field can be represented as a linear combination 

n 

ajO~j, aj e Z 2 . (68) 

Then, we solve the equation (J53jl c\ ^ = x f° r the n basis elements c Ktfl , 

k = ai, a n assigning in an arbitrary way the signs ±1 to the square root a/ X (/« 2 a0- 
This means that there exist 2 n different sets of {c K ^}, and thus 2 n different operators 
(for a fixed value of /i) , which can be precisely parameterized as in Eq. f)3U|) . Once the 
signs of c K ^, k — (Ti, a n , are fixed, the rest of the 2 n — n coefficients c K ^ can be found 
using the expansion (jo^j) and the relation (J2llj) in the form c K+a ^ = c K ^c a ^x {~^olk) 
leading to the following result: 

C-a,[j, C a i<7i+..+a n o- n ,/i C aia - lj ^ t C a2 (T2+..+a n ov l ,/iX (^1^1 ^fll^l ■■ Q"n,0~n) AO 

(n— 1 n \ 

fc=i j=fc+i / 

and c 0)/J = 1. 

To illustrate how this procedure works let us apply it to the case of GF{2 2 ). We 
choose the normal, self-dual, basis {9, 9 2 } (see Appendix A) in GF(2 2 ), where 9 is a 
root of the primitive polynomial x 2 + x + 1 = 0, so that 9 3 = 9 + 9 2 . The solution of 
Eq. (jHH) for, say /i = 9 3 = 1, is (below we will omit the index /i in the coefficients c QjAt ) 

c = 1, c e = ±i, c 6 2 = ±i. 

Then, the last coefficient is given by 

cg3 = c e+9 2 = c e c 2x (-9) = (±i) (±i) (-1) , 

and one can see that there exist 4 different possible operators Vg3. A similar calculus 
can be made for the operators Vq and Vg2. 

It is convenient to fix positive signs of the coefficients c KjAt corresponding to the 
elements of the field basis, i.e. c a;(TiiA1 = a/x ( a ?°f A 4 )) I = 1,-,ti, and form the "first" set 
of the rotation operators, with these coefficients . Then, all the other sets of V^ v can 
be obtained according to ()3U|). 

Once we have the coefficients c QjAt one can easily obtain the corresponding phase 
factors for the displacement operator (J50|) . For instance, fixing the "first" set of rotation 
operators in the above example as 

Vq = diag(l, 1, i, — i), Vp. = diag(l, i, 1, — i), Vg3 = diag(l, i, i, — 1), (69) 

we obtain the following phase factors appearing in the displacement operator: 

<i>(9,9) = i, <KM 2 ) = i, <l>(o,e 3 ) = i, 

<P(9 2 ,9) = 1, <P(9 2 ,9 2 )=z, <P(9 2 ,9 3 )=z, 
(f,(9 3 , 9) = -i, cf)(9 3 , 9 2 ) = -i, <P(9 3 , 9 3 ) = -1. 
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Another set of rotation operators can be obtained, for instance, keeping Vg and Vgi as 
in the above and multiplying the operator Vg3 in by Xg\ 

Vg — diag(l, 1, i, — i), Vg2 — diag(l,i, 1, — i), Vgs — diag(l, — i, i, 1), 

which leads to some changes in the phases of the displacement operator: 

0(0,0) = -*, 0(M 2 ) = i, 0(0,0 3 ) = *, 

0(0 2 ,0) = 1, 0(0 2 ,0 2 )=*, 0(0 2 ,0 3 )=*, 
0(03, 0) = -i, 0(0 3 , 2 ) = -i, 0(0 3 , 03) = L 

10. Appendix C: Determination of / (/^, //) 

In practice, the function / (/i, //) is determined by solving the following equation 

c a ,f,c a ^ = x (af (ji, A*')) c ^+m' ( 70 ) 

for all values of a belonging to the basis of the field. The rest of the elements of the 
field do not provide additional information due to (j2HJ). 

As an example we consider the case of GF(2 2 ). Following the general procedure, 
we choose the self-dual basis (0, 2 ) in the field and fix the rotation operators as in (fMJ) . 
The equation (|7U|) for / (0, 0) has the form 

c a ,gc a ,g = x (af (0, 0)) c a , , 
and for different values of the parameter a we get 

« = 0, 1-1 = X (0/(0,0))1, 
a = 2 , i-i = x(9 2 f(9,0))l, 

leading to the only possible solution / (0, 0) = 2 . 
Similarly we obtain for the rest of the / (//, //): 
for/(0,0 2 ) 

« = 0, i-t = x {of{e,e 2 ))t, 

a = 6\ *-l = X (0 2 /(0,0 2 ))*, 
so that / (0, 2 ) = 0, and following the same idea we can determine 

/ (0, 3 ) = 2 , / (0 2 , 2 ) = 0, / (0 2 , 3 ) = 0, / (0 3 , 3 ) = 3 . 
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11. Appendix D: Another symplectic operators 

11.1. U operator 

In order to introduce another operator with similar properties to V^, let us define 
that two lines are orthogonal if the states corresponding to these ones are related via 
the Fourier transform 

Then, there exists an operator U^k) = FV^k) such that its geometrical application 
rotates the line corresponding to the state \k) (conjugate to and then, transforms 
the rotated line into an orthogonal one. This operator is obtained from as 

U ll = FV^ = C -*M( K \> ( 71 ) 

K£GF(d) 

so that, U^ZaU^ = Z a for all a, p. The action of U p on the operator X@ can be obtained 
using the operational relation in (fTTj) and the relation ((7j) 

U^XpUl = F^ZpV^F^ = exp(-^ (/?, p))Z^X^ (72) 

which means that the ^-transformation also represents a sort of rotation: operators 
transform eigenstates of the set of displacement operators labelled with points of the 

ray [3 = pa 

{-^j Z ai X pai , Z^X^^ , . .. } (^3) 

to the eigenstates of the set labelled with points of the ray (3 = {p + p') 1 a 

{I, Z( 1+wt /) ai X MQ!1 , Z(i + ^i) a2 X^ a2) ...} . (74) 

Similarly, as we cannot reach the ray a = using operators, one cannot reach the 
ray (3 = using operators. Having both operators U^i and we can transform any 
ray into any other. 

Let us find the Wigner function of a state transformed by and U p operators. 
For fields of odd characteristic, p ^ 2, one can use the explicit form ffTTTTf) of and after 
some algebra we obtain 

W p -(a,(3) = W p (a,[3-pa), p = V ll pV^. (75) 

In the same manner we evaluate the Wigner function of a state transformed by U p 
operator, p ^ 2: 

W~ p {a,(3) = W p {a-p(3,(3), ~p = U llP Ul (76) 

In other words, the transformation of a state by V p and/or operators leads to a 
covariant transformation of the Wigner function for the fields of odd characteristics. 
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11.2. S operator 

The last operator which transforms the operators labelled with points of one line into the 
operators labelled with the points of some other line is the so-called squeezing operator, 
which in the basis of eigenstates of Z a operators has the following form 

K£GF(p n ) 

It is easy to see that 

W- p (a,P) = W p (Za,C 1 P), P = S( P S$. 

The squeezing operator naturally appears as a result of consecutive application of V p 
and U u : 

w P > = w SipsV 

where /i is an arbitrary element of GF (p n )*, p ^ 2. 
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